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We present preliminary results for a high statistics study of the topological charge distribution in 
the SU (3) Yang-Mills theory obtained by using the definition of the charge suggested by Neu- 
berger fermions. We find statistical evidence for deviations from a gaussian distribution. The 
large statistics required has been obtained by using PCs of the INFN-GRID. 
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1. Introduction 

We perform a high statistics study of the topological charge distribution in the SU{'i) Yang- 
Mills theory adopting the definition suggested by Neuberger fermions as it is discussed in a series 
of recent studies [l]-[4]. The numerical studies have been initiated in ref.[8, 9] and completed in 
ref. [10] where a systematic study at different volumes and values of the lattice spacing has been 
performed in order to obtain a precise and reliable determination of the topological susceptibility. 
Their distributions are obtained by a population between 1500 and 3000 topological charge for each 
volume. This already impressive number of configurations allowed them to measure the variance 
of the distribution which is given {Q^) and the topological susceptibility at a ~ 5% level but it 
was not enough to show non-gaussianity. Moreover in [10] it has been numerically confirmed that 
a physical volume larger than (Ifm)"* guarantees that finite-volume effects are negligible at their 
statistical precision. 

The aim of this study is to look for non-gaussianity in the topological charge distribution of the 
SU{3) Yang-Mills theory. Here we study three lattices at the same physical volume ~ (1.12fm)^ 
with about ten time more statistics than in ref. [10] in order to emphasize the deviations from the 
gaussian distribution. We stress that in order to search for such very small subleading effects it is 
necessary to be sure that all the systematics of the calculation cannot either simulate or hide the 
effect and therefore the only reliable theoretical framework is the one provided by the topological 
charge definition suggested from Ginsparg- Wilson fermions. 

This challenging Monte Carlo calculation has been recently made possible by a few important 
improvements in the algorithmic sector which guarantee the reliability and the feasibility of high 
statistics [1 1]. Moreover algorithms for zero mode counting with no contamination from quasi zero 
modes, optimized to run fast on a single processor, are now available. The great computer effort 
has been performed in the frame of the INFN GRID project which allowed us to use the computer 
resources shared in the scientific Italian network provided by INFN along this year. 

2. Theoretical procedure 

In the following we use the standard plaquette action of the SU (3) gauge field. The massless 
lattice Dirac operator D satisfies the Ginsparg-Wilson relation 

Y5D + DY5=dDY5D . 

In this frame the topological charge density can be defined 

q{x) = -^Tr[Y5D{x,x)] (2.1) 

where d = a/{l -\- s), a is the lattice spacing and the shift parameter s has been fixed in our calcu- 
lation at the value s = 0.4. The topological charge is obtained from the lattice by computing for 
each configuration the difference between the numbers of zero modes with positive and negative 
chiralities that is the index v of the Dirac operator: 

V = — ?i_ , 
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which is directly related to the topological charge Q: 

V = Q = a'^Y^q{x) 

X 

As it is stressed in ref.[8] the asymptotic distribution of the topological charge follows a gaus- 
sian form centered in zero, with the only parameter {Q^) as variance : 

^ ,exp(-;^). (2.2) 



Corrections to this formula are expected to be suppressed by terms of 1 /Nc'^ and 1 /V where V is 
the lattice volume and Nc is the number of colours of the theory. 



3. Our project 

In the Table 1 it is shown the summary of the features of our three runs altogether with the 
information about the occupancy of the Grid network. The number of configuration aie written 
in form of a product of the number of simultaneous runs sent to the Grid network and the total 
number of configurations planned for each run. The fourth column indicates the computer time 
taken for each configuration running on a single processor. The fifth column contains the number 
of configuration sequentially processed in a single run sent to the Grid network. In the case of the 
largest volume the desired number of configurations (30000) was not achieved at the time of the 
Conference. 



Lattice 


/3 


ro/a 


Total # Confs 


Time 


# Confs for each run 


Ram (MB) 


124 


6 


5.368 


58x600 


1 h/conf 


10 


256 


144 


6.0938 


6.263 


120x250 


4 h/conf 


4 


512 


184 


6.2623 


8.052 


200x150 


12 h/conf 


2 


1100 



Table 1: Summary of our lattices, physical volume ~ (1.12/m)4, ro — 0.5fm. The run of the 184 lattice was 
not completed at the time of the presentation. 

We stress that all our runs are performed over a single processor using its local memory. The 
parameters of the lattices are chosen to fix the physical volume at the value ~ (1.12/m)4 at three 
lattice spacing in order to estimate the size of the discretization effects. 

In order to check the statistical independence of the measurements we have first studied the 
topological charge autocorrelation function. The worst case of the 1 84 lattice is shown in Fig. 1 . 
This figure shows a typical autocorrelation time of order of about 40 in units of Monte Carlo 
updates. Each update is composed of one heat-bath and 6, 7, 9 cycles of over-relaxation of all 
the links for each volume. We have chosen to separate the subsequent measurements by 500, 
1500 and 2500 updates for the volumes 124, ^44 ^§4 respectively. This large separation among 
measurements does not cost too much in computer time and allows us to consider the measurements 
as statistically independent. 
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Figure 1: Autocorrelation of the topological charge as function of computer updates in the worst case of the 
largest volume IS'*. 

4. Distributions at large statistics: preliminary data 





Data from ref. [10] 


Our preliminary data 




L/a 


# Conf s 






# Conf s 


6 


12 


2452 


1.63(8) 


1.638(13) 


34800 


6.0938 


14 


1405 


1.54(6) 


1.565(13) 


30000 


6.2623 


18 






1.442(18) 


13561 



Table 2: Comparison with ref. [10] for the first two lattices. The three physical sizes have about the same 
value 1.12 fm. 

In the Table 2 we show our preliminary data for the correlator {Q^) in comparison with the 
results of ref. [10]. The values are in complete agreement each other and moreover it is clear that 
the errors scale with the 1 / \/N law as it is expected due to their statistical nature. 

In order to give a graphic representation of our data, following [10], we show in Fig. 2 the 
histograms of the topological charge of the first two lattices for which the runs were completed at 
the time of the presentation. We underline that the topological charge is a discrete variable and 
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therefore the gaussian form must not be confused with the normal distribution even if the mathe- 
matical expression is the same. Here plotting the histograms of the topological charge values in 
comparison with the continuos line of the gaussian, we want only to give a graphical representation 
of the data and the comparison with the asymptotic formula must be done comparing the value at 
the center of the column with the level of the column. The border at the top indicates the statistical 
error (1 sigma) of the column. At first glance it is worth to note that the data at Q = are sensitively 
higher than the theoretical gaussian value. 




Figure 2: Histograms of the topological charge for 12^ and 14^ lattices. The border at the top indicates the 
statistical error of the column. In each picture the value of the is reported. 

To be quantitative we test the hypothesis that large statistic data are compatible with the lead- 
ing asymptotic distribution eq.(2.2) taking as parameter the measured value of {Q^). We use the 
test of value that is reported in each of the pictures in Fig. 2. The value of for the data at 
large statistics is greater than the critical value value ^ corresponding to more than five standard de- 
viations therefore the data are not compatible with the pure gaussian distribution of the topological 
charge. 

In order to study the behaviour as function of the number of collected topological charges we show 
in Fig. 3 for the smallest volume 12^ the histograms at increasing statistics. We plotted two his- 
tograms of 1000 and 10000 values. The test indicates that the first sample is compatible with 
a pure gaussian within one sigma while in the second sample the data show a larger x^ value cor- 

^The critical value for a distribution witli 12 degrees of freedom is 14, 30, 52 for 1, 3, 5 standard deviations 
respectively. 
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responding to a deviation of more than three sigmas and therefore the test-hypothesis should be 
rejected. 




Figure 3: Histograms with 1000 and 10000 configurations for the 12^ run. In the first sample the data are 
statistically compatible with a pure gaussian while for the second one the data hegins to show a deviation of 
more than 3 sigmas from a pure gaussian. 

The most interesting result we have to show is the value of the fourth cumulant of the topolog- 
ical charge distribution: 

This quantity is zero for a pure gaussian distribution and the previous data from ref. [10] are con- 
sistent with zero. The values we found are definitely different from zero. We quote preliminarly: 

K4(latticel2^) = 0.609 ± 0.093 and K4(latticel44) = 0.442 ± 0.078 

where the enors are 1 jacknife standard deviation. 

5. Outlook 

We plan to finish the runs of the third volume in a few weeks reaching the prefixed value 
of 30000 configurations, then we will complete the checks in particular the dependence from the 
volume size, for which we plan to run a test volume at 16^, j8 = 6 with L ~ 1.5fm. 

Even though the analysis presented and the results are still preliminary they already shown 
that the topological charge distribution presents deviations from gaussianity. This goal has been 
reached with a statistics of 30000 configurations. 
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